How weak values emerge in joint measurements on cloned quantum systems 
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A statistical analysis of optimal universal cloning shows that it is possible to identify an ideal 
(but non-positive) copying process that faithfully maps all properties of the original Hilbert space 
onto two separate quantum systems. The joint probabilities for non-commuting measurements on 
separate clones then correspond to the real parts of the complex joint probabilities observed in 
weak measurements on a single system, where the measurements on the two clones replace the 
corresponding sequence of weak measurement and post-selection. The imaginary parts of weak 
measurement statics can be obtained by replacing the cloning process with a partial swap operation. 
A controlled-swap operation combines both processes, making the complete weak measurement 
statistics accessible as a well-defined contribution to the joint probabilities of fully resolved projective 
measurements on the two output systems. 
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Quantum mechanics cannot be explained by a simple 
realistic model of physical properties, because it is impos- 
sible to jointly determine the values of non-commuting 
observables. If a final measurement on an initial state 
| ip) determines the values of an observable A with eigen- 
states | a), it is impossible to know what the result for 
an observable B with eigenstates | b) would have been 
for the same system, and vice versa. Specifically, a mea- 
surement of B performed after the measurement of A 
will not have the same result as a measurement of B per- 
formed before the measurement of A, since the projection 
of the quantum state onto an eigenstate of A completely 
changes the statistics of B. 

In measurement theory, the projection of the state is 
related to decoherence in the interaction with the meter 
system. One way to avoid the changes in the statistics 
of B is therefore to weaken the measurement interac- 
tion until the decoherence effects are negligible. This ap- 
proach to quantum measurement was introduced in 1988 
by Aharonov, Albert and Vaidman, and is now widely 
applied to the study of quantum systems [l - 11 1 . In weak 
measurements, the value of an observable A between the 
preparation of an initial state | ip) and a final measure- 
ment of | b) is obtained from the small shift of a meter 
caused by a weak interaction. The average result ob- 
tained from many measurements of | ip) with the same 
final outcome | b) is given by the weak value, 
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The right hand side of this equation shows the relation 
between the complex weak value and the eigenvalues A a . 
The analogy with conventional quantum statistics sug- 
gests that the weak value is an average of the eigenvalues 
A a for a complex conditional probability of the eigen- 
state | a) in the pre- and post-selected ensemble defined 
by | ip) and | b). This interpretation is also consistent 



with the notion that probabilities should be defined as 
expectation values of the projection operators | a) (a | , 
since the complex conditional probabilities in Eq.(JT|) are 
given by the weak values of | a) (a |. Significantly, the 
negative real parts of these weak conditional probabili- 
ties can explain quantum paradoxes, as demonstrated in 
a number of recent experiments [2h10||. 

Based on the statistical interpretation of weak values, 
it is possible to define a joint probability of the measure- 
ment outcomes | a) and | b) [I3 - H6j |. Specifically, the 
joint probability of a and b is obtained by multiplying 
the conditional probabilities of a with the appropriate 
probability of b given by | (6 | -0) | 2 to obtain 
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Weak values therefore suggest that the correct joint prob- 
ability of a and b is obtained by using the expectation 
value of the operator product of the projectors on | a) and 
| b) . Remarkably, this definition of complex joint prob- 
abilities was already introduced by Kirkwood in 1933 as 
an alternative to the Wigner function, and later gener- 
alized to discrete systems (l7Hl9j |. For arbitrary mixed 
states, the complex joint probability p wea ^(a, b) provides 
a complete characterization of the density matrix for any 
pair of basis sets {| a}} and || b)} with non-zero over- 
laps between all | a) and | b) [12, EH, Since weak 
values thus correspond to a fundamental expression of 
quantum statistics that was already studied long before 
weak measurements were introduced, one might expect 
that other implementations of joint measurements should 
reveal the same correlations between non-commuting ob- 
servables. In this letter, it is shown that the joint proba- 
bilities derived from the statistical interpretation of weak 
measurements also appear in the correlations between 
cloned quantum systems, indicating that cloning can be 
used to confirm the statistical correlations between non- 
commuting observables described by weak values in a 
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physical context that is completely different from weak 
measurements. 

Although cloning usually focuses on the attempt to 
copy the quantum state, the statistical interpretation of 
quantum mechanics implies that this can only be done by 
copying the actual physical properties, whether defined 
by the state or not. Ideal universal cloning should there- 
fore copy all physical properties faithfully to both out- 
puts, resulting in perfect correlations between the clones 
for all measurements. Of course, such an ideal cloning 
process is impossible, as stated by the well-known no- 
cloning theorem [20| . However, approximate cloning can 
be optimized by exploiting the symmetry of the quantum 
states [21| . In the following, it is shown that the errors in 
the optimal cloning process can be interpreted as an ac- 
cidental distribution of white noise to one of the outputs, 
leaving a contribution that defines ideal cloning in terms 
of a non-positive linear map. Significantly, this ideal 
cloning process describes perfect correlations between all 
physical properties of the output systems, supporting the 
notion that the cloning process actually copies all phys- 
ical properties equally, whether they are eigenvalues of 
the input state or not. 

Since the physical properties of the output clones in the 
ideal cloning term are perfectly correlated, separate mea- 
surements performed on the two clones can be interpreted 
as simultaneous measurements on the same system. The 
joint probabilities obtained from the clones therefore rep- 
resent the correlations between non-commuting proper- 
ties in the single input system before the cloning process. 
Significantly, the analysis of the ideal cloning term repro- 
duces the joint probabilities previously observed in weak 
measurements, as given by Eq. ([2]). The only signifi- 
cant difference is that cloning does not define a temporal 
sequence, resulting in equal contributions from the se- 
quence (a,b) and (b,a). Since these contributions have 
opposite imaginary parts, only the real part of the com- 
plex joint probability P we ak( a > ^) appears in the corre- 
lations between the clones. The imaginary parts can be 
obtained by a partial swap operation between the input 
state and white noise, which differs from the cloning pro- 
cess only in the phase of the superposition between swap 
and no swap. It is therefore possible to recover the tem- 
poral order of weak measurements from the direction of 
the partial swap operation that can be used to realize the 
optimal cloning process. 

The starting point of the following analysis is the for- 
mulation of optimal quantum cloning as a map between 
the single d-dimensional Hilbert space of the input and 
the bosonic (or positive parity) states of the two output 
systems given by (| ip;m) + \ m;ip)) [io|, [2l|. By mixing 
over a complete orthogonal basis set { | to) } , all informa- 
tion other than the input state | ip) is eliminated. The 
output density matrix of the optimal cloning process can 
then be written in terms of the contributions from | ip) in 
output 1, from | ip) in output 2, and from the coherence 



between the two, 
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Here, {| to)} can be any orthogonal basis of the single 
system Hilbert space. The first two terms can be inter- 
preted as valid quantum operations in their own right, 
corresponding to the random distribution of the input 
state and a completely random state to the outputs with 
a total probability of d/(d + 1). The remaining output 
probability of 1) is then provided by the coherence 

term, which describes a non-positive map of the input to 
the two output systems. 

In its most compact form, the contribution of the co- 
herence term to the output can be represented by an 
operator C^, that is not self-adjoint and hence describes 
complex-valued contributions to the statistics, 
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By itself, this operator defines a non-positive linear map 
of the input state | ip) (ip | that faithfully copies the input 
state to both outputs. Specifically, the local states ob- 
tained when tracing out one of the two output systems 
are both given by | ip), 
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Thus the non-positive state has a perfect cloning fi- 
delity of one. Moreover, this state also describes perfect 
correlations between all quantum fluctuations of | ip). 
If the same projective measurement onto an orthogonal 
basis {| a)} is performed in both systems, only con- 
tributes to the probabilities of finding the same result in 
both outputs, as shown by the projection 
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This is the signature of a physical copy of the input sys- 
tem: even an unknown property a is faithfully duplicated. 
It may therefore be justified and useful to think of the 
quantum coherence in the optimal cloning process as 
a duplication of all physical properties in Hilbert space, 
where the quantum coherences in the initial Hilbert space 
arc transferred to directly observable correlations be- 
tween the separate Hilbert spaces of two systems. 

If different measurements are performed in the two out- 
put systems, the contributions of to the joint proba- 
bilities of | a) and | b) are given by 



pc(a, b) = (a; b \ | a; b) 

= {iP\b)(b\a){a\iP). 
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As explained in the introduction, this is exactly equal 
to the complex joint probability P wea k( a '^) m Eq.©, 
which is obtained from a weak measurement of | a) (a \ 
followed by a final measurement of | b)(b \. Specifically, 
the weak values of | a) (a | followed by a post-selection of 
| b) correspond to the conditional probabilities obtained 
by dividing pc(&, b) by the total probability of b given by 
|(6 | "0)| 2 - The weak value of A can then be expressed as 
a conditional average of the joint probability pc(a, b), 
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The statistical interpretation of weak values is therefore 
consistent with the correlations independently obtained 
from quantum cloning. 

It might seem a bit strange that the ideal cloning pro- 
cess G,p corresponds to a specific sequence of a and b. In 
the real cloning process given by Eq.(J3]), there is no pre- 
ferred sequence since the output includes an equal mix- 
ture of the ideal cloning processes and correspond- 
ing to the sequences (a, b) and (&, a), respectively. As a 
consequence, the imaginary parts of the joint probabili- 
ties pc(a,b) cancel, and only the real parts show up in 
the experimentally accessible joint probabilities obtained 
from an optimal cloning process. Specifically, the actual 
output probabilities p(a, b) measured in the output of the 
optimal cloning process are given by 

p(a,b) = (a;b\E clone (\^)(^\)\a;b) 

= (|<a | i>)\ 2 + \(b | ^)\ 2 + 2Re (p c (a, &))) .(9) 

In addition to the real part of the complex joint prob- 
ability pc(a,b), cloning errors result in a "background" 
obtained from the probabilities |(o | ?/>)| 2 and \{b | ip)\ 2 - 
Intuitively, this corresponds to the random results for a or 
b obtained when the cloning process distributed the white 
noise input to one of the outputs. Since the marginal 
probabilities p(a) and p(b) obtained for each system can 
be used to determine \{a | ip}\ 2 and \(b | ij))\ 2 , it is a 
straightforward matter to subtract the background noise 
and to determine the real part of pc{.a>, b) from p(a, b). 

For completeness, it would also be desirable to obtain 
the imaginary part of the cloning statistics. In weak 
measurements, the change from real to imaginary weak 
values is achieved by replacing the weak measurement 
operator with a weak unitary transformation generated 
by the projector on | a) 0, [ll IH-El • Effectively, this 
corresponds to a phase change in the superposition be- 
tween no measurement and projective measurement that 
characterizes the weak measurement interaction. In the 
present case, a similar effect can be achieved by changing 
the phase in the superposition between the distribution 
of the initial state to system I and the distribution to 
system 2. In terms of quantum logic operations, the dis- 
tribution of the input can be achieved by performing or 



not performing a swap operation. The eigenstates of this 
operation are the states with positive (bosonic) and neg- 
ative (fermionic) parity under the swap operation. The 
imaginary superposition of swap and no swap results in 
a unitary operation with eigenvalues of (I + i)/^/2 and 
(I — i)/V% for the positive and negative parity states. 
The operation thus corresponds to a partial swap repre- 
sented by the square root of the unitary operator for the 
complete swap. The positive linear map of this partial 
swap or root-swap operation on an input product state 
of | i])) and the maximally mixed state I/d is given by 
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Thus the imaginary part of the expectation values of 
given by pcifl,b) represent the (real) correlations that 
build up between the two systems during the continuous 
transfer of quantum information from one system to the 
other in a unitary swap operation. Essentially, the cor- 
relations of the partial swap operations provide a com- 
plete map of the dynamical structure of Hilbert space 
that complements the statistical structure observed in 
the correlations between clones (23^ . 

Finally, it might be worth noting that an efficient com- 
bination of optimal cloning and partial swaps can be real- 
ized by a controlled-swap operation with a control qubit 
input in a superposition state of | Z }+ | l z ). Any super- 
position of swap and no swap can then be accessed by 
simply measuring the corresponding superposition in the 
control qubit output. Since is a direct representation 
of the coherence between swap and no swap, its contri- 
butions to the joint probabilities of a and b can be ob- 
tained from the differences between the joint probabilities 
obtained for opposite control qubit coherences. Specifi- 
cally, the real part of pc{ a , b) can be obtained from the 
difference in the joint probabilities for the positive and 
negative superpositions of | Z ) and | l z ) given by the 
A-basis eigenstates | 0^) and | l^), while the imaginary 
part is given by the differences between the correspond- 
ing y-basis eigenstates | Oj,) and | 1^). The quantum 
controlled-swap then permits complete quantum tomog- 
raphy of the input state based on two non-commuting 
von Neumann measurements performed on the separate 
output systems and two different measurements on the 
control qubit output, 

Re(pc(a, b)) = d (p(0 Xl a, b) - p{l x ,a, b)) 
lm(p c {a, b)) = d (p(0 y ,a, b) - p(l y , a, b)) . (II) 

The complete density matrix of an unknown input state 
p can then be reconstructed directly by 
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The quantum controlled-swap thus maps the complete 
quantum coherence of the input onto the correlations 
between a and b in the output. In contrast to weak 
measurements, no approximate limit is required and the 
statistical structure of Hilbert space appears as a well- 
defined part of conventional measurement statistics. In 
the controlled-swap, the subtraction of background noise 
is achieved by taking the difference between the proba- 
bilities of opposite measurement outcomes for the control 
qubit, so complex probabilities are observed directly as 
polarizations of the qubit in the ATF-plane of the Bloch 
sphere. Note that this procedure is technically similar 
to the observation of the complex wave function in the 
polarization of the meter photon reported in How- 
ever, the controlled-swap achieves much higher visibili- 
ties, since the measurements are not performed in the 
weak interaction limit. 

The results presented above indicate that the com- 
plex joint probabilities obtained by multiplying the non- 
commuting measurement operators provide a consistent 
representation of statistical correlations between non- 
commuting measurements in quantum mechanics. In 
particular, the non-positive linear map represented by 
appears to be the quantum mechanical analog of a 
faithful copy of the input system to both output systems. 
Since this linear map works for all input states and does 
not depend on the choice of measurements performed on 
the output, the origin of the correlations between the 
measurement outcomes can be traced back to the single 
Hilbert space of the input state. Thus, quantum cloning 
and partial swaps provide an approach to joint measure- 
ments that avoids many of the ambiguities of sequential 
or joint measurements performed on the same system. 

It is especially remarkable that the results correspond 
to those observed in weak measurement, since the physics 
of the measurement procedure are obviously quite differ- 
ent. While weak measurements distinguish between the 
weakly measured result a and the post-selected result b, 
cloning and partial swaps are completely symmetric in 
the output measurements. The measurement interaction 
is well separated from the cloning process, so the correla- 
tions observed cannot be induced by the dynamics of the 
measurement interaction. The observation of weak mea- 
surement statistics in cloning and in partial swaps there- 
fore supports the interpretation of weak measurements 
as back-action free measurements of joint probabilities 



performed jointly. The statistics obtained from the joint 
measurements correspond to those obtained in weak mea- 
surements and thereby support the interpretation of weak 
measurements in terms of joint and conditional probabili- 
ties. A statistical analysis of quantum mechanics in terms 
of complex joint probabilities and non-positive maps can 
thus help to uncover unifying principles in quantum in- 
formation and improve our understanding of the physics 
behind non-classical effects. 
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In conclusion, the effects of quantum coherence in 
optimal cloning processes can be described by a non- 
positive linear map that represents the quantum mechan- 
ical analog of a perfect copying process. This perfect 
copying process maps the statistical correlations of non- 
commuting measurement in a single Hilbert space onto 
two separate systems, where both measurements can be 
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